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DEVELOPMENT OF A SIMPLIFIED PROCEDURE FOR ROCKET 
ENGINE THRUST CHAMBER LIFE PREDICTION WITH CREEP 

SUMMARY 


An analytical method for predicting engine thrust chamber life 
is developed. The method accounts for high pressure differentials 
and time-dependent creep effects both of which are significant in limiting 
the useful life of the shuttle main engine thrust chamber. The hot-gas- 
wall ligaments connecting adjacent cooling channel ribs and separating 
the coolant flow from the combustion gas are subjected to a high pressure 
induced primary stress superimposed on an alternating cyclic thermal 
strain field. The pressure load combined with strain-controlled cycling 
produces creep ratcheting and consequent bulging and thinning of these 
ligaments. This mechanism of creep-enhanced ratcheting is analyzed 
for determining the hot-gas-wall deformation and accumulated strain. 
Results are confirmed by inelastic finite element analysis. Fatigue 
and creep rupture damage as well as plastic tensile instability are 
evaluated as potential failure modes. It is demonstrated for the NARloy Z 
cases analyzed that when pressure differentials across the ligament 
are high, creep rupture damage is often the primary failure mode for 
the cycle times considered. 



DEVELOPMENT OF A SIMPLIFIED PROCEDURE FOR ROCKET ENGINE 
THRUST CHAMBER LIFE PREDICTION WITH CREEP 

1.0 INTRODUCTION 


Life predictions of regeneratively liquid-cooled rocket engine thrust 
chambers have usually been based on low cycle fatigue evaluations. 
Experiments, however, have shown that low-cycle fatigue is not always the 
controlling failure mode. The chamber coolant walls progressively thin 
and bulge with each firing cycle until failure occurs due to tensile 
instability. In a recent report for NASA [1]* an analysis of the 
deformation and strain accumulation due to plastic ratcheting of the hot- 
gas-wall due to successive firings of a cylindrical test chamber was 
performed. The tensile instability mode of failure was analyzed as well 
as low-cycle fatigue. A simplfied analytical procedure for predicting 
thrust chamber life was developed. 

The analyses of [1], however, were for an experimental cylindrical 
thrust chamber where the pressure differential between the coolant channel 
and the combustion gas chamber is significantly lower than the corres- 
ponding pressure differential in the space shuttle main engine (SSME). 

In addition, the experimental engine firing period is limited to several 
seconds, whereas the SSME engine firing time is measured in minutes. 

The procedure developed in [1], therefore, does not include any creep 
and is applicable only for short cycle times. For the SSME thrust chamber 
where cycle times are much longer, the time-dependent response of the 
material must be included in the creep ratcheting analysis. The higher 
pressure differential in the SSME thrust chamber must also be considered. 
The present work is aimed at developing a simplified method for predicting 
the number of times the SSME can be fired before the liner ligaments fail. 
The method would allow the designer to optimize the various engine design 
parameters to achieve longer life. The accuracy of the simplified method 

*Numbers in brackets refer to References at end of Report. 
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is verified via the finite element technique by performing one or two 
cycles of inelastic analysis on a simple ligament beam model. 

The scope of the work is as follows: 

1. - Development of the creep-ratcheting solution for determination 
of ligament (hot-gas-wall) deformation and strain accumulation. 

2. - Analysis of failure modes and criteria. Develop failure 
criteria by investigating: 

(a) fatigue 

(b) creep rupture damage, and 

(c) plastic/creep tensile instability 

3. - Synthesize a simple design procedure for predicting thrust 
chamber life based on 1 and 2 above and provide sample calculations. 

Each of these items will be discussed in the sections to follow. 

The assumptions made in the development of the analysis procedure for 
each item are discussed in their respective sections. 

2.0 DEVELOPMENT OF THE CREEP RATCHETING SOLUTIONS FOR DETERMINATION OF 
HOT-GAS-WALL DEFORMATION AND STRAIN ACCUMULATION 

2.1 Description of the thrust chamber 

The experimental plug nozzle thrust chamber is shown in Figure 1. 

This figure shows a typical plug nozzle assembly consisting of the 
contoured centerbody and flanged cylinder along with the cross-sectional 
details of the cylinder. The inner liner of the cylinder contains axial 
flow coolant channels similar to those of the SSME thrust chamber. The 
inner wall separating the coolant flow and the combustion gas is subjected 
to severe thermal cycling during startup and shutdown. 
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The ligaments in the SSME thrust chamber are subjected to higher 
differential pressure between the coolant channels and the combustion 
gas side of the liner and longer firing cycles than the experimental 
chamber. The high primary stress due to this differential pressure is 
superimposed on the alterlating cyclic thermal straining. The pressure 
load combined with strain-controlled cycling produces creep ratcheting 
and consequent bulging and thinning of the ligaments. This mechanism of 
creep-enhanced ratcheting is analyzed to develop an analytical model for 
determining the hot-gas-wall deformation and accumulated strain. 

2.2 Inelastic Strains 


The hystersis loop for a typical loading cycle, not including the 
initial chilldown, is shown schematically in Figure 2. OA and AB represent 
the initial elastic and plastic portions during heat up. Creep occurs 
during operation at elevated temperature causing the stresses to redistribute. 
The stress field generated after severe plastic straining is compressive 
at the peak of the transient and tends to relax asymptomatically to the 
stationary state related to the sustained primary load due to pressure. 

The relaxation of the thermally induced stresses due to creep results in 
the creep-enhanced ratcheting mechanism mentioned above. The creep relaxation 
is depicted by the vertical line BC in the figure. CD and DE represent 
the elastic and plastic portions during cool down. Due to creep relaxation, 
initial yielding occurs earlier resulting in a larger plastic portion. 

Thus, as a result of including creep, the elastic portion is reduced and 
the plastic portion increased. 


Denoting the average temperatures of the ligament and closeout wall 
by T. and T Q , respectively. The inelastic strain range in the hoop direction 
due to differential thermal expansion EB and creep relaxation BC, is given 
by: 


Ae 


Pi 


i a i 





/E+ Aa c /E 


( 1 ) 
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FIGURE 2 TYPICAL LOADING CYCLE LOOP 
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where a. and a Q are the thermal expansion .coefficients of the ligament 

and closeout, respectively, and (T^ - T Q a 0 ) max and (T^- T 0 a 0 ) m -j n are 

the maximum and minimum thermal strains, respectively, that occur during 

the loading cycle. A a is the relaxation stress and S and S are 

c c ^max y min 

respectively the ligament material absolute yield strengths corresponding 

to the ligament average temperatures at the time in the cycle when 

(T a - T a ) and (T. a . - T a ) . are calculated. The inelastic 
i i o o max l 1 o o min . 

compressive strain generated within the EBC part of the cycle is reversed 

by plastic straining along the CDE line. 

Though Ac 1 as given by Equation (1) accounts for the major portion 
Pi 

of the inelastic strain range in the hoop direction, there is also a tempera- 
ture drop across the ligament which causes bending, since the ligament 
ends are constrained. This thermally induced bending, though acting only 
for a short portion of the cycle, may enhance the ratchet strain. Its effect 
may be assessed by computing the elastic energy of the thermally induced 
bending stresses and correcting the hoop strain accordingly. Conservatively 
assuming that all of the available elastic energy goes into plastic straining 
of the ligament, the correction due to thermally induced bending is given 
by [1]: 


Ac 


ii 

P 


E(aAT) 2 

i " T2('l-v) 2 S y 


( 2 ) 


where AT is the temperature drop across the ligament. 

The total inelastic strain in the hoop direction for the complete 
cycle is thus: 

e ' ■ 2 K, * *1) (3) 
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2.3 Creep relaxation 


For low pressures, as in the case of the experimental thrust chambers, 
the relaxation of stress can be determined by using the uniaxial solution 
for zero total strain rate. Assuming Norton's law governs creep flow; 

e = Ba r ( 4 ) 

where B and r are material constants, the governing differential equation 
is : 


1 da , D r „ 
idt + Bo =0 


(5) 


On integrating, the stress after relaxation within time t is given 


by 


a = a p [l + T(r-l)] 1/1 - r 


( 6 ) 


where a 

T 


= peak stress, and 
r- 1 

= BEOp t denotes nondimensional 


time 


Here, a = S^, the yield stress. Thus the creep relaxation stress is given 
by: 



(7) 


Equation (6) has been verified by finite element analysis and found to 
be very accurate for predicting the mid-plane relaxation stresses for the 
low pressure cases. When the pressures are high, however, the relaxation 
proceeds faster and reaches the steady state much more rapidly. Figure 3 
shows the mid-plane stress relaxation curves for the case of no pressure 
(which is coincident with the finite element results for the low pressure 
case of 547 psi) as given by Equation (6) and the case of high pressure 
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FIGURE 3a STRESS RELAXATION CURVES (MID-SURFACE) 



(p = 3000 psi) as obtained from finite element analysis described later 
in the report. It is seen that for the high pressure case the mid-plane 
stress relaxes to almost zero and hence the creep relaxation stress is 
given by : 


Aa c = S y (8) 

For intermediate pressure cases, the relaxation stress may be obtained 
by interpolation from Figure 3. For conservative evaluation. Equation 
(8) may be conveniently used in all pressure cases. 


The relaxation curves plotted in Figure 3a have been obtained by titnrig 
the creep law constants for NarloyZ at 1325°F. Stress relaxation curves 
that are applicable in the general case can be generated from Equation (6) 
for the zero pressure condition by plotting the dimensionless stress para- 
meter a/a vs. nondimensional time r for various r. These curves are 
shown in Figure 3b and may be used for estimating Ao c when the pressure 
differentials are lower than 1000 psi. For higher pressure differentials, 
Eq. (8) is applicable in the general case. 
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2.4 Incremental deformation of ligament 

The corresponding curvature and hence the deflection is determined 
using the method developed in [1]. The method uses a yield surface for 
combined bending and membrane loading to determine the incremental inward 
bulging and progressive thinning near the center of the ligaments at the 
inner liner of the thrust chamber. The derivation details are given in 
[1] and hence only an outline is presented here. 

2,4.1 Mechanics of inelastic distortion 


The approach used to determine the strain increments within the load 
cycles is explained in Figure 4. The beam of Figure 4a simulates the 
response of the ligament subjected to pressure-induced bending and cyclic 
straining in the hoop direction. The yield surface for a beam subjected 
to bending and hoop force is given by [2]: 


m + n 2 - 1 = 0 


(9) 


where m and n are dimensionless variables defined by: 


m 


n 


M 

M ( 

N 


(10a) 

(10b) 


In the above, N and M denote the hoop force and bending moment, while N Q 
and M q denote the yield hoop force and yield bending moment given by: 


N o = 2H S y 


M rt = H 2 S 
o y 


(Ha) 

(Hb) 


for a rectangular beam of unit width, height 2H and yield stress S^. 

If shear is also included in the model, the yield surface is three- 
dimensional as shown in Figure 5 and is given by [3]: 
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X - Hoop Direction 
Y - Axial Direction 
Z - Radial Direction 


b) 

Plastic Contract 
tion of Liga- 



Bending Component 6 

Plastic Stretching 
of Ligament 

Known hoop 
component X 
of resulting 
plastic strain 
vector 

o -N/(2H S y ) 
m = M/M q “ M/ (H* S y ) 

X - 2 Hi 

§ » u*i 


FIGURE 4 PLASTIC RATCHET IN A BAR OF RECTANGULAR 
CROSS SECTION SUBJECTED TO SUSTAINED 
BENDING MOMENT AND CYCLIC HOOP STRAIN 
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FIGURE 5 YIELD SURFACE FOR COMBINED BENDING, SHEAR AND MEMBRANE LOADING 


- k = 0 


( 12 ) 


where 


k =vl - s : 


2t 

and s = -F- is the dimensionless shear stress. 

y 

Using the normality law, the generalized strain rates are: 


i = Af = 2An 


(14a) 


i.xfE.AVTr? 


(14b) 


^ = = A 2s - 


1 - s : 


(14c) 


where ft is an arbitrary positive scalar and A, e, are the generalized 
hoop, curvature and shear strain rates respectively. The relationship 
between the generalized strain rates and the hoop strain rate, curvature 
rate and shear strain rate are, respectively [1]: 


A = 2He! 


(15a) 


e = h 2 k 


(15b) 


<P = Hy 


(15c) 


From Equations (3), (14) and (15), the curvature and shear strain are, 
respectively: 
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(1- S 2 £i 


Y= &■£ ? 


\ >FV \" ) 

The curvature and shear strain as determined from Equation (16) and 
(17) for each cross-section can then be integrated along the length of 
the ligament to obtain the corresponding bending deflection 6 i and shear 
deflection 62 • 

2.4.2 Creep-induced deflection 

In addition to & 1 and $ 2 , the steady creep strain due to pressure 
also contributes to the total deflection. This additional deflection 
is determined by treating the ligament as a clamped beam under uniform 
loading, as described below. 

Consider the beam of rectangular cross-section shown in Figure 6 . 

Then Bernoulli's 1 approximation of plane sections gives for the strain rate 
at point P in the x-direction: 

S ' - a^t ' '«■ (,8 > 


Using Norton's Creep Law: 


i = Ba r 
x x 


where B and r are material constants, gives: 


1 ru\ 1 /' r 1 / r 

°x ^77 (v > 2 

The bending moment M is given by: 


M = / Z o x dA 
A 
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SECTION THROUGH P 


FIGURE 6 CO-ORDINATE SYSTEM AND BEAM DIMENSIONS 
FOR DERIVATION OF EQUATION (22) 
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where A is the cross-sectional area of the beam. Substituting from (20), 
the following differential equation is obtained for the deflection of the 
beam: 


= » r (22 > 
Following [4], for the ligament under pressure loading. Figure 7: 



where m is a constant to be determined from the boundary conditions. Using 
(23) and the clamped boundary conditions i.e., v'= 0 at x =0, l , gives 
the maximum deflection: 


v (r) = 
max v ' 



where F(r) has been tabulated for various r in Table I. 
the derivation of (24), the reader is referred to [4]. 


(24) 

For details of 


Thus total deflection per cycle 

6 = 61 + 62+63 (25) 

2.4.3 Ligament distortion 


Experimental evidence shows that the deformed shape of the ligament 
can be approximated by a linear variation in thickness as shown in Figure 8 . 
Thus using the linear thinning model developed in [1], the thinning after 
N cycles is then given by: 


_ N 6 w 
N (£ + w) 


(26) 


where £ is the width of the ligament in the hoop direction and w is the 
width of the rib. 
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Uniformly loaded, clamped beam co-ordinate system. 


FIGURE 7 LIGAMENT UNDER PRESSURE LOADING 
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3.0 ANALYSIS OF FAILURE MODES AND CRITERIA 


The ligaments are subjected to incremental permanent deformations 
during each firing cycle of the thrust chamber. The geometry of these 
ligaments changes as the incremental strains accumulate. They are subjected 
to incremental bulging and progressive thinning near the center of the 
ligament. Fatigue and creep rupture damage as well as plastic tensile 
instability are considered as potential failure modes and evaluated 
herein. 


3.1 Fatigue 


The strain range at the minimum ligament section increases with 
progressing distortion and thinning. For fatigue calculations, the maximum 
local hoop strain range at the minimum ligament section can be obtained 
from the average hoop strain range by integration if the geometry of the 
distorted ligament is known. The method of determining the maximum local 
strain range in the minimum ligament section is the same as that developed 
in [1] and hence not repeated here. For the linear thinning model of 
Figure 8b, the effective strain range in the minimum ligament section for 
entering the fatigue curve is then: 


+ Ei • E2 • + e 2 • 

^•imn z min 2 min 


where £ i m j n > the local hoop strain in the minimum ligament section is: 



and q is the exponent in the Ramberg-Osgood type stress strain relation 
e = Aa q . The average hoop strain is given by: 

£l avg = a(T i " V (29a ) 

The axial strain E 2 min in the minimum ligament section is also given by: 
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(29b) 


where is the average temperature of ligament and T Q is the average 
temperature of closeout. 

For the linear variation assumed, the thickness t m _. and t max after 

h 

the N cycle are given by: 


_ 2H(& + w) - N 6 w 
min U + w) 


(30) 


= 2HU + w) 2 + N 6 ft w 
r max (X+wJ 7 


(31) 


where 5 is the deformation per cycle from (25). 

Equations (27) through (31) along with the fatigue curve are then 
used to determine the fatigue life. Note that though the procedure parallels 
that in [1], the resulting strain range is not the same, since the deflection 
and hence the thinning per cycle is now larger due to creep effects. 


3.2 Creep Rupture Damage 

The compressive stress field generated at the peak of the transient 
relaxes during the creep period to the steady state pressure stress. When 
the pressure acting on the ligament is low, as in the case of the experi- 
mental thrust chamber, the finite element results show that there is negligible 
redistribution of the bending stress due to pressure and the hoop stress 
generated by the thermal transient remains compressive up to the end of 
the creep period. Thus, creep rupture damage is not a potential failure 
mode for such cases. 

When the pressure acting on the ligament is higher, as in the SSME 
thrust chamber, the finite element results show that the hoop stress relaxes 
quite rapidly becoming tensile on the hot-gas-wall surface of the ligament. 

The pressure stress is redistributed during the creep period and relaxation 
occurs to this steady state value. Thus creep rupture damage must be 
included in the failure analysis for high pressure cases. 
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The stress relaxation curves for the hot-gas-wall surface of the ligament 
are shown in Figure 9 for the two pressure cases discussed. Interpolation 
between these curves may be used to obtain the relaxation curves for 
intermediate pressure cases. Also, the curves can provide an estimate 
of the time over which the stress remains tensile for consideration in 
the creep rupture damage analysis. 


The creep rupture damage evaluation is based on the tensile stress 
at the end of the creep period. This may or may not be the steady state 
stress. Since the creep period is known, the time for the pressure stress 
to redistribute must be less than this creep period for it to reach its 
steady state value. This time may be estimated by using Calladine's work, 
[5]. He has derived a simple formula for the time taken for the redistri- 
bution of stress in a structure from an initial state to a final steady 
state to be practically complete. The time taken is expressed in terms 
of the time taken for creep strain to become equal to a certain multiple 
of the elastic strain in a constant stress creep test performed at a stress 
level which depends on the steady state stress in the structure. This 
time for redistribution is given by: 


t 


ss 




(32) 


where 

°ss = steac *y state stress 

e = elastic strain 

B and r are constants in the creep law. 

Though the above formula was derived for kinematically determinate structures, 
it provides a reasonable estimate for the redistribution of mid-channel 
pressure stress in the present case. 


The finite element results for the clamped beam model show that for 
the creep exponent r = 6.65, the mid-channel steady state stress at the 
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FIGURE 9 STRESS RELAXATION CURVES (HOT-GAS-WALL SURFACE) 
FOR NARLOY Z FROM FINITE ELEMENT ANALYSIS 
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surface is approximately 0.8 of the initial bending stress. Calladine in 
[6] has studied several structures subject to creep for determining the 
way in which the greatest stress in the structure varies with the creep 
exponent r when the geometry and the loading remain unchanged. He showed 
that a linear interpolation between the cases r = 1 (elastic) and r = » 
(perfectly plastic) may be used for estimating the greatest steady state 
stress for arbitrary r. Following this thought, the steady state stress 
in the present case may then be estimated by: 


°ss ' Fa in 


(33) 


where a. is the initial stress and F is plotted for various r in Figure 10. 
This plot has been generated by connecting the two known values of F at 
r = 1 and r = 6.65. The maximum initial elastic stress can be determined 
using, for example, finite element solution. Its value can also be 
estimated with sufficient accuracy by treating the ligament as a clamped 
beam under pressure loading. 


Having determined the steady state stress from (33), this is substituted 
in (32) to estimate redistribution time t . If t is less than the cycle 
creep period, creep rupture damage evaluation is based on o ss - If, however, 
t is greater than the cycle creep period, creep rupture damage is usually 
small and may often be excluded from failure evaluation. Since the stress 
on the hot-gas-wall surface must become tensile before reaching the steady 
state value, a case could conceivably arise where the stress has become tensile 
but did not quite reach a ss by the end of the creep period. In order to 
provide conservative evaluation that would preclude such cases, it is 
possible to determine a fraction of t which should be compared to the 
cycle time period. Based on the analysis performed, the curves of Figure 9 
indicate that for cycle times limited to 0.1 t , contribution of creep 
rupture damage can be conservatively ignored in the failure analysis. 


26 




The creep rupture damage is determined from the stress to rupture 
curves. Since the thickness of the ligament changes with each cycle, a 
numerical procedure is used to evaluate the damage. Neglecting the effect 
of moment redistribution along the ligament, the initial pressure induced 
stress for each cycle is based on t - n from Equation (30). The total 
fatigue and creep rupture damage can be determined by linear summation, 
strain range partitioning or other damage summation methods. 

Kachanov's creep rupture damage theory is described in Appendix A. 

If the stress rupture curve (log a vs, log t ) is not a straight line for 
a given temperature then damage accumulation will depend on the stress 
level and there could be a large deviation from unity if Robinson's life 
fraction rule is used in a linear manner [7]. Kachanov's theory can be 
useful in explaining such deviations. In the present case, either theory 
can be conveniently used for evaluating creep rupture damage. 


3.3 Plastic Instability 

At the peak of the heat transient, the ligament yields in compression. 
Thus, the stresses are compressive during the creep period. The ligament 
then yields in tension at the cold end of the cycle and ductile failure 
due to plastic tensile instability must be considered. The critical 
effective strain and critical ligament thickness were developed in [1] 
and are given by: 


— _ 2q(l-a+a 2 ) 2 

°cr " (2-a) 


(34) 


t 


cr 


2He“ < * 


(35) 


where a = -^ and ai and o 2 are the principal stresses in the hoop and axial 
directions, respectively. 
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4.0 THRUST CHAMBER LIFE PREDICTIONS 


The analysis method developed herein provides the cyclic strains 
needed to perform a fatigue damage evaluation and the stresses needed to 
evaluate creep rupture damage. The latter is quite significant for the 
high pressures used in the SSME thrust chamber. The resulting combined 
damage may limit the life of the thrust chamber, i.e., the number of times 
it can be fired before ligament cracks can be anticipated. This analysis 
includes the effect of progressive geometry changes of the ligament on 
strain ranges and maximum stresses. 

Plastic tensile instability is an independent mode of failure which is 
simultaneously checked with progressive deformation as an alternative pos- 
sibility of failure. This mode of failure appears to be limiting for 
certain materials such as OFHC copper. NarloyZ appears to exhibit cyclic 
hardening mechanisms which when combined with corresponding shift of yield 
frame related to kinematic hardening results in gradually diminishing thin- 
ning in consecutive cycles. The interaction between creep hardening 
mechanisms and cyclic hardening is, of course, quite complex and beyond 
the scope of this evaluation. The present analysis does not take account 
for these complexities of material response and therefore may predict 
premature plastic tensile instability where such hardening is significant. 
Fatigue and creep rupture damage may, in fact, be more limiting than 
plastic tensile instability for NarloyZ. 

The results obtained using the analyses described herein can be used 
for conservative estimates of the life of the thrust chamber. If the 
number of cycles to failure based on fatigue and creep rupture damage of 
the gradually thinned ligament are determined to be less than those 
needed to reduce the ligament thickness below the critical value resulting 
in plastic instability, then the failure mode is obvious and the cycles to 
failure due to fatigue and creep rupture damage correspond to the 
estimated life of the thrust chamber. 
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If, however, the number of cycles to reach the critical thickness 
obtained from the tensile instability analysis is determined to be less 
than that due to fatigue and creep rupture damage, then although the 
evaluation remains conservative the dominating failure mode for materials 
such as NarloyZ is not obvious. In these cases, plastic instability 
analysis and creep and fatigue damage evaluation provide lower and upper 
bounds on the cyclic life of the thrust chamber, respectively. Such situa- 
tions may occur when the pressure differential across the ligament is low, 
as for example in test chambers but they are basically not expected for 
the SSME thrust chamber where the pressure differentials are substantially 
higher. 

Numerical calculations for various cases of NarloyZ using the 
procedure developed herein are presented in Appendix B. 
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5.0 FINITE ELEMENT ANALYSES 


A finite element model of the ligament was formulated to confirm the 
analytical procedure developed herein. The cross-hatched area in Figure 
11 shows the ligament portion modeled. The finite element mesh along with 
the co-ordinate system and boundary conditions is shown in Figure 12. Two- 
dimensional isoparametric solid elements were used for the model. Assuming 
generalized plane strain conditions, the analysis was performed for NarloyZ 
using the thermal history obtained from [1] and shown here in Figure 13. 

The material was assumed to be elastic perfectly plastic and creep was 
allowed to occur for four minutes at the hot end of the cycle. Norton's 
law was used for the creep flow. The pressure load was taken as 3000 psi 
(20.68 MPa) and acted over the ligament portion shown in Figure 12. 

The finite element mesh with the node numbers is shown in Figure 14. 

The results of the finite element analysis are shown in Figures 15 through 
18. The stresses are plotted for nodes 133 to 143 which lie along the center- 
line of the ligament in the thickness or Y-direction. Sx is the hoop stress 
and SIGE is equivalent stress. Figure 15 shows the stresses due to a pressure 
of 3.000 psi (20.68 MPa). The thermal transient is then applied and the 
ligament yields compression, as illustrated in Figure 16. During creep, the 
hoop stress relaxes quite rapidly becoming tensile at the inside surface in 
about one second, as shown in Figure 17. Relaxation occurs to the steady 
state stress as shown in Figure 18. 

The results for the high pressure case show that the stress is redis- 
tributed during creep and the steady state stress at the hot-gas-wall surface 
is approximately 80% of the initial value. Moreover, since the stress 
on the hot-gas-wall surface becomes tensile very rapidly compared to the cycle 
time, creep rupture damage is a potential failure mode that should be 
included in the failure analysis. 
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FIGURE 16 PEAK OF THERMAL TRANSIENT 

(Y =0 is Hot-gas-wall Surface) 
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To investigate if such is also the case with low pressures, a finite 
element analysis was performed for NARloy-Z using a pressure of 547 psi. 

(3.77 MPa). The finite element mesh, the model boundary conditions and 
the thermal history were the same as that for the high pressure case. 

The results of this analysis are shown in Figures 19 through 21. The 
stresses are again plotted along the centerline of the ligament in the 
thickness or Y-direction. Figure 19 shows the stresses due to pressure 
only. Figure 20 shows the stresses at the peak of the transient. The ligament 
has yielded in compression. Figures 21 and 22 show the stresses after creep 
for one second and 240 seconds respectively. It is seen that there is 
negligible stress redistribution and the hoop stress remains compressive 
even after four minutes as compared to the previous case of high pressure 
(p = 3000 psi (20.68 MPa)) for which the hoop stress became tensile, at the 
inside surface after only about one second. 

The analysis for the low pressure case was then extended to include 
creep for eight minutes. It was found that the hoop stress remained com- 
pressive, indicating that creep rupture damage does not occur for low pressure 
cases . 

The stress relaxation curves for the middle and inside ligament surface 
are plotted in Figures 3 and 9, respectively, for the two pressure cases 
analyzed. 

The above analyses indicate that creep effects are significant for 
the SSME. thrust chamber, and creep rupture damage should be considered 
since pressures are high. For the experimental thrust chamber, however, 
creep has little effect on the cycles to failure. 
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FIGURE 19 INITIAL CONDITION - P = 547 PSI ONLY 
(Y =0 is Hot-gas-wall Surface) 
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FIGURE 21 CREEP FOR ONE SECOND DURING THERMAL TRANSIENT - AT = 865°F 
(Y =0 is Hot-gas-wall Surface) 





FIGURE 22 CREEP FOR 240 SECONDS DURING THERMAL TRANSIENT - AT = 865°F 
(Y = 0 is Hot-gas-wall Surface) 




6.0 COMPARISON OF RESULTS 

The numerical results obtained using the analytical procedure developed 
herein are compared with the finite element results in Table II. The test 
case analyzed was the experimental thrust chamber geometry with a pressure 
differential of 3000 psi (20.68 MPa) corresponding to the SSME thrust chamber. 
The two are in very good agreement indicating that the analytical method is 
modeling the physical behavior, of the ligament accurately. 

TABLE II 


COMPARISON OF ANALYTICAL VS. FINITE ELEMENTS RESULTS 
Residual Deflection After One Cycle 


Analytical 

Results 

Finite Element 
Results 

Deflection of Pres- 
sure Surface/Cycle 

Deflection of Pres- 

| sure Surface/Cycle 

! 

.34 mils (.0086mm) 

.34 mils (.0086 mm) 


7.0 CONCLUSIONS 

An analytical method for predicting engine thrust chamber life is 
developed herein. The method accounts for high pressure differentials 
and time- dependent creep effects both of which are significant in limiting 
the useful life of the SSME thrust chamber. 

It is demonstrated from the cases analyzed that when pressure differ- 
entials across the ligament are high, creep rupture damage is often the 
primary failure mode for the cycle times considered. It is therefore 
important that the thrust chamber is designed to minimize the combined 
stress and temperature effects on progressing creep damage of the material. 

The simplified method of analysis developed herein provides the strain 
ranges and stresses needed to perform fatigue damage and creep rupture damage 
eval uations. 
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The simplifed method of analysis developed herein provides the number 
of firings which will result in a plastic tensile instability failure but 
takes no credit for the plastic and creep hardening mechanisms in materials 
such as NarloyZ. These hardening mechanisms have not been quantified* 
and may greatly extend the number of firings which may be tolerated prior 
to tensile instability. 

The simplified method of analysis developed herein provides an excellent 
tool for making a parametric study to determine the coolant channel design 
and operating parameters which would maximize the number of repeated firings 
that could be achieved prior to failure. 


*The relevant material properties could be measured on simple specimens 
but complex temperature and strain histories would have to be used in order 
to obtain the properties needed for analysis purposes. 
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APPENDIX A 


KACHANOV'S BRITTLE CREEP RUPTURE DAMAGE THEORY 


The usual approach to evaluating creep damage accumulation is based 
on Robinson's life fraction rule: 

(A-l) 

where t^ is the time at a given stress level and t pi - is the rupture time 
for that stress. The concept implied in the model of Equation (A-l) is 
that damage accumulation is a linear process independent of the order 
of the succession of applied stress. 



An alternative approach to evaluating brittle creep rupture damage 
is that due to Kachanov. His theory relates the damage rate to the maximum 
tensile stress in the current undamaged net area of the ligament cross- 
section. If A denotes the total cross-sectional area and A r denotes 
the currently effective load carrying area, the creep rupture damage can 
be defined as: 


D = 



(A-2) 


Kachanov postulated the following power relation for the creep damage 
rate: 


dD 

it 


= C(a r ) 


v 


a \V 




(A-3) 


where a and a r are the initial and current stresses corresponding to A 
and A r , respectively, and C and v are material constants. D = 0 corresponds 
to material in the virgin state and D = 1 corresponds to failure. 
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For a constant stress level. Equation (A-3) gives on integration: 

' - ‘V+V ■ Co v t + * (A— 4 ) 

where K is the constant of integration. 

Using the conditions that when t = 0, D = 0 and when t = t , D = 1 gives: 


I = 



(A-5) 


and 


r v - 1 

t r (v+l) (A-6) 

Thus Equation (A-4) becomes 


1 



If the plot of log o vs log t r for a given temperature is a straight line, 
then v in Equation (A-7) is independent of the stress level and can be 
obtained from the slope of this line. In this case a single damage accumu- 
lation curve is obtained as shown schematically in Figure A-l and the 
failure criterion is identical to the life fraction rule even though the 
damage accumulation process is nonlinear. 

If the above condition does not hold for a particular temperature, then 
damage accumulation will depend on the stress level which in turn governs v. 
When v is stress-sensitive, there could be a large deviation from unity if 
the time ratio t/t r is used in a linear manner, as illustrated in Figures 
A-2(a) and A-2(b). Kachanov's theory can be potentially useful in explaining 
these deviations in such situations. 


In the present case, however, either theory can be conveniently used 
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FIGURE A-l DAMAGE VS. TIME RATIO WHEN v 
IS INDEPENDENT OF STRESS LEVEL 



OAM AGE. D 


(a) 



Vi < v 2 
cr(vi) < a(v 2 ) 


Program: a 2 oi 

Result: £ (|-) > 1 
i = 1 r 1 



FIGURE A-2 DAMAGE VS. TIME RATIO WHEN v IS 
DEPENDENT ON STRESS LEVEL 




to evaluate the creep rupture damage, the two being related through 
Equation (A-7). Figure A-3 is a schematic illustration of D representing 
the area damage using Kachanov's theory and t/t representing the linear 
damage theory vs. cycle ratio. This figure graphically illustrates the 
physical difference between these two theories. 
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DAMAGE, D 
TIME RATIO t/t 



0 1.0 


CYCLE RATIO N/N Q 


FIGURE A-3 SCHEMATIC VARIATION OF D and t/t r VS CYCLE RATIO 
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APPENDIX B 
NUMERICAL EXAMPLES 


Example I ; 

Consider the ligament geometry shown in Figure B-l. The material 
is NARloy Z, with the following properties: 

Yield Stress S = 24,500 psi 

Young's Modulus E = 16.35 x 10 6 psi 
Poisson's Ratio v = 0.34 

Coefficient of Thermal Expansion, a = 9.5 x 10 -6 in/in°F 


For a unit width of the ligament Ap = 3000 Ib/in. 

From Equation (l), taking a. = a = a and fs + S )/2 = S , gives the 

u vy max ^min 7 y 

inelastic strain range due to differential thermal expansion and creep 
relaxation: 



a 




- J o) 


min 



(B-l) 


From thermal analysis: 


(T. - T ) - fl 

r. - t . Y 

\ l o/max \ 

i o min/ 


865 °F 


Creep law: 


e = Ba r 

The ligament temperature was taken to be 1325°F, with the following creep 
constants estimated from Table B-l: 

B = 5.5225 X TO" 7 
r = 6.65 

and a is in ksi. 
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a = 
2H = 
w = 


0.0664" (1.686 mm) 
0.035" (0.889 mm) 
0.05" (1.27 mm) 


FIGURE B-l LIGAMENT GEOMETRY 












The mid-plane stress relaxes to almost zero as seen from Figure 3. Thus 
from Equation (8) 


A a = Sy = 24,500 psi 


From Equation (B-l) 



= 9.5 x 10 


-6 


x 865 - 


2 x 24500 1 
16.35 x 10 b J 


L 24500 
. "16.35 x 10 b 


= 0.006719 


The temperature drop across the ligament aT ~ 210°F. Thus from Equation (2) 
the correction to the plastic strain range due to thermally induced bending: 

_ 16.35 x 10 6 x (9.5 x 10" 6 x 210 1 2 
Pi 12(T~- 0.34 ) 7 x 24500 


= .000508 

From Equation (3), total inelastic strain 

ci = 2 ( .006719 + .000508) 

= .014454 


From Equation (16), taking n ~ k gives: 


K = 


- 1 _ £i 


R " JT 


where the radius R is assumed for the length shown in Figure B-2 
Thus 

r = .0175 1 pi in 

.014454 L ' dL1 
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Deflection 61 = 2R(1 - cose) 


= 2 ( 1 .211 )[1 - {( 1 . 211) 2 - (.0166) 2 }*/1.211] 
= .0002276" 


Equation (17) is approximated by: 


Y = 2sci 


where the dimensionless shear stress s is given by [ 1 ] : 


Thus, s 


s = -Rsn (! - x ) 

3000 / .0664 

.0175 x 24500 \ 2 


= .2323 - 6.9971x 


(B-2) 


Substituting in Equation (B-2) integrating and determining the shear 
deflection at x = H/2 due to pressure loading, gives 


6 2 = .000112 


From Table I, for r = 6.65, F(r) = .1168 by linear interpolation 


Thus from Equation (24) for cycle time t = .0667 hr (4 min) and pressure 
of 3 ksi, the creep-induced deflection: 


6 3 


5.5225 x IQ - 7 x .0667 x (. 06641 2 (l x .0664 2 \ 6 ' 65 
4 x .0175 \l 6 x .61757 


.1168 


= . 0000002 " 

Thus total 6 = 61 + 62 + 6 3 
= .0003398" 
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From Equation (26), thinning after N cycles: 


t ,. 0003398 x .05N 
L N (.0664 + .05) 

t = .000146 N inches ( B_3 ) 

N 


Plastic Instability : 

The strain hardening parameter q is approximately given by LeRc as 
'S - S \0.6 


= 0.2 


For S = 55 ksi and S = 24.5 ksi 
o y 


. _ A o i 55 - 24.5 
q - °‘ 2 1-243 


0.6 


= .228 


From Equation (35) 


t = ( .035 )e 
cr 


-0.228 


= .02786 


Thinning for instability = .035 - .02786 

= .00714" 


( B-4) 


Equating (B-3) and (B-4) 

N = 49 cycles. 

Fatigue and Creep Rupture Damage 

The procedure of section 3.0 was used for evaluating fatigue and 
creep rupture damage. The fatigue curve of Figure B-3 and the stress to 
rupture curve of Figure B-4 were used. as the basis of this evaluation. 
Linear summation was employed to estimate total damage. 
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From Figure 10, for r = 6.65, factor F = 0.8. Thus from 
Equation (33) the steady state stress: 


a 


ss 


O.S o jn 


( B— 5 ) 


The initial stress is estimated by treating the ligament as a clamped 
beam under pressure loading. From Figure B-l the mid-surface length of the 
beam = 0.0914". For simplicity and conservatism, the loading is assumed 
to be over this engire length.* Thus 


D- 


in 


3000 ( . 09 14 ) 2 
4~f . 035 ) 2 


= 5115 psi 


From Equation (B-5) a ss = 0.8 x 5115 = 4092 psi 


Substituting in Equation (32) 

t = 1 „ 2.3 5.115 

ss 5.225 x 10“' x (4 .092 ) b - 5 x 6.65 x T6.35 x 10* 

= . 0167 hrs. = 60 seconds 

The creep constants for 1325°F were used in the finite element analysis 
and hence are also used here in evaluating t $$ . Since 0.1 t gs is less than 
the creep period of four minutes, creep rupture damage is based on a . 

The FORTRAN program listed in Appendix C was used for performing the 
necessary calculations. The input and output obtained are enclosed in 
this Appendix. 


The notation for the listed results is as follows: 
NCYC - Number of cycles 

TMIN - t mi - n , minimum ligament thickness - in. 


♦Assuming a .0914" long beam that is loaded over a partial length of 
& = 0664" gives a slightly lower stress. 
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TMAX - t , maximum ligament thickness - in. 
max 

EEQVT - effective strain range % 

N/Np - fatigue damage 

STRESS- primary stress due to pressure 

T/T D - creep rupture damage 

TOTAL USEAGE - cumulative usage factor 

NF CYCLES - cycles to failure 


The output shows that since T/Tp is substantially higher than N/Np, 
creep rupture damage is the dominating failure mode, fatigue being negligible. 
The estimated life is determined to be nine cycles. 
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CONSTANTS input 


Q« 


0,22800 


DELTA 


0.33980E-03 

INCHES 

E1AVG 


0*82175 

PERCENT 

H 


0.1750OE-01 

INCHES 

L 


0.66400E-01 

INCHES 

W 


0.50000E-01 

j 

INCHES 

P 


3000.0 

LBS/INCH 

T/CYC 


4.*0000 

minutes 

SF 


0,80000 



FATIGUE CURVE 


CYCLES 

10.000 

J00.no 

aooo.o 

0.5000oE*06 


STRAIN RANGE 

18.000 

5.9000 

\.0000 

0.10000E+00 


STRESS^TD-RUPTURE curve 


HOURS 

o^ioonoE+oo 

O.BOOftO 

8.0000 

100.00 

tooo.’o 


STRESS (KSI) 

6.0000 

4.0000 

2.5000 

1.5000 

1.0000 



RESULTS FOP EaCh CYCLE 


NCYC 

TM f N 

T M a y 

FEOVT 

N/NO 

STRESS 

T/TO 

TOTAL IJsAoE 

CYCLES 

\ 

o,3R8S4E-oi 

0.3S063E-0I 

1.6554 

0.71265E-03 

9.1261 

0.97710E-01 

0.98922E-01 

10 

2 

0,397PflE-01 

0.3S167E-01 

1 .6675 

0.72352E-03 

9.1609 

0.10201 

0.20115 

9 

3 

0 r 39562E-01 

n.35250E-0t 

1 .6798 

O.73964E-03 

9.1961 

0.10651 

0.30890 

9 

a 

0 r 3fl«lbE-0l 

0.35333E-QI 

1 .6922 

0.74603E-03 

9.2317 

0.11129 

0.92039 

9 

5 

0,3927oE-01 

0 # 35 A 1 6E-Q 1 

t .7099 

0.75769E-03 

9.2679 

0. 11620 

0.53739 

a 

6 

0*3<I1?«E-01 

0.3S50ftE-0l 

1.717ft 

0.76963E-03 

9.T045 

0.12190 

0.65951 

6 

7 

0,3397bE-01 

0 . 35563E-Q l 

1.7309 

0.78187E-03 

9.3915 

0. 12665 

0.78719 

b 

8 

0 ^3383?E-0 1 

0.35666E-01 

1.7492 

0.79440E-03 

9.3791 

0.13258 

0.92051 

6 

9 

0.33bflfcE-01 

0.357<|9E-01 

1.7577 

0.80723E-03 

9.9171 

0.13859 

1 .0599 
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Example II 


Consider the ligament geometry shown in Figure B-6. which corresponds 
to that of the SSME thrust chamber. The material is NARloy Z, with the 
following properties: 

Yield Stress = 24,500 psi 

Young's Modulus E = 16.35 x 10 6 psi 
Poisson's Ratio v =0.34 

Coefficient of Thermal Expansion, a = 9.5 x 10 -6 in/in°F 


For a unit width of the ligament Ap = 3000 lb/in. 

From Equation (l), taking a. = a = a and + S )/2 = S , gives: 

vy max ^min' y 



OL 


(r 


i 



r o) 


min 



(B-l ) 


From thermal analysis, for 120% power level, 


(t. - T ) - (' 

r. - T • ) 

\ l o/max \ 

i o min/ 


1500°F 


Creep law: 


e = Bo r 

For an average ligament temperature of 1000°F, the creep constants from 
Table B-l are: 

B = 6.2 x 10 -16 
r = 12.5 

and o is in ksi. I 
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£ = 0.04" (1 -016 mm) 

2H = 0.028" (0.711 mm) 

w = 0.045" (1.143 mm) 


FIGURE B-6 LIGAMENT GEOMETRY 




The mid-plane stress relaxes to almost zero. Thus from Equation (8) 


A o = Sy = 24,500 psi 


From Equation (B-l) 



9.5 x 10" 


x 1500 - 


2 x 24500 1 
16.35 x 10 b . 


, 24500 

16.35 x 10 b 


= 0.0127515 


The temperature drop across the ligament aT - 650 0 F * Thus from Equation (2): 

Ar" = 16.35 x 10 6 x (9.5 x 10~ 6 x 650)2 
Pi 12(1 - 0.34) z x 24500 


= 0.0048681 


From Equation (3), total inelastic strain 

d = 2( .0127515 -+ .0048681) 

= .0352391 

From Equation (16), taking n ~ k gives: 


K = — = — 
K R H 


where the radius R is assumed for the length shown in Figure B-2 
Thus 

R = = .03973" 


.0352391 
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Deflection = 2R(1 - cose) 


= 2 ( . 3973 ) C 1 - {(.3973 ) 2 - (.01) 2 } */.3973] 
= .0002517" 


Equation (17) is approximated by: 


where 
Thus, s 


Y = 2 s ex 


S = Tf ^“ (f ' x ) 


HS > 

3000 


.014 x 24500 


- 


.1749 - 8.7464X 


( B-2 ) 


Substituting in Equation (B-2) integrating and determining the deflection 
at x = 1/2 gives 


62 = .0001233 


From Table I, for r = 6.65, F(r) = .1168 


Thus from Equation (24) for cycle time t ■= .0667 hr. 


63 


= 6.2 x 10 - 16 x .0667 x (.04) 2 / 3 x .04 2 \ 6,65 
4 x .014 yl 6 x .OUJ 


.1168 


» 0 

Thus total £ = 61 + 5 ? + 63 
= .000375" 
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From Equation (26), thinning after N cycles: 


t = .000375 x .045N 
N (.04 + .045) 

t^ = .0001985N inches 


Plastic Instability: 


For S = 55 ksi and S = 24.5 ksi 
u y 



From Equation (35) 


t = (.028)e“ 0 - 228 
cr 


= .0222915 


Thinning for instability = .028 - .0222915 

= .0057085 


(B-4) 


Equating (B-3) and (B-4) 

N = 29 cycles. 

Fatigue and Creep Rupture Damage 

The procedure of section 3.0 was used for evaluating fatigue and 
creep rupture damage. The fatigue curve of Figure B-3 and the stress to 
rupture curve of Figure B-4 were used as the basis of this evaluation. 
Linear summation was employed to estimate total damage. 
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From Figure 10, for r = 12.5, factor F = 0.77. Thus from Equation (33) 
the steady .state stress: 


The initial stress is Estimated by treating the ligament as a clamped 
beam under pressure loading. The mid-surface length of the beam = 
0.0625". For simplicity and conservatism, the loading is assumed to 
be over this entire length. Thus 


o 


in 


3000 (.0625) 2 
4 ( .028 ) 2 


= 3737 psi 


From Equation (B-5) a s$ = 0.77 x 3737 a 2878 psi 


The ligament in the present example has a large temperature variation 
across its thickness. The surface temperature of the ligament on the 
hot gas side is approximately 1325°F while that on the coolant side is 
about 675°F. If the creep constants are based on an average temperature 
of 1000°F in determining t , then from Equation (32) 

t = 1 x 2 .3 x 3.737 

SS 6.2 x 10 _1G x (2.878) 12 - 5 12.5 16.35 x 10 3 


= 4.46 x 10 8 seconds 


In this case 0.1 t is greater than the cycle time of four minutes 
and the stress remains compressive at the end of the creep period. Thus, 
there is no creep rupture damage and the life estimate is based on fatigue. 
The FORTRAN program listed in Appendix C was used for performing the necessary 
calculations. The input and output are included in this Appendix. Since 
thinning cannot continue beyond t , the fatigue life is estimated to 
be 206 cycles. 
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If on the other hand, a more conservative approach is employed and 
the creep constants for the surface temperature of 1325° F are used in 
determining t $$ (since creep rupture damage occurs at the surface) then 
from Equation (33): 

a = .8 a. = 2990 psi 
ss in p 

From Equation (32): 

t 1 x 2.3 x 3.737 

SS 5.225 x 10" 7 x (2.99) 6 ' 65 6.65 16.35 x 10 3 

= .0984 hrs. = 354 seconds 

In this case 0.1 t is less than the cycle time of four minutes. 

The stress at the surface thus becomes tensile and creep rupture damage 
must be included in the life evaluation. The necessary calculations were 
again performed using the FORTRAN program listed in Appendix C. The input 
and output are included in this Appendix. 

The output shows that creep rupture damage is now the dominating 
failure mode, fatigue being relatively small. The estimated life is only 
twenty cycles. 

This example illustrates the sensitivity of the results to the creep 
properties. It also helps explain the wide discrepancies in the failure 
life of thrust chambers observed during tests, which have been attributed 
to local hot spots. The example here shows that if some ligaments do get 
hotter than others, creep rupture damage could cause early failure. In 
other instances if the ligaments remain cooler, the chamber lasts much 
longer. 
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cunstants input 


Q« 


0,22600 


DELTA 

a 

0.37500E-03 

INCHES 

F J A VG 

a 

1.‘4250 

PERCENT 

H 

a 

0, 14000E-01 

INCHES 

L 

a 

0.40000E-01 

INCHES 

W 

a 

0.45000E-01 

INCHES 

P 

a 

0.0000 

LBS/INCH 

T/CYC 

a 

0.0000 

MINUTES 

SF 

a 

o.'oooo 



fatigue curve 


CYCLES 

10,000 

100.00 

4000,0 

0,'50000E*0b 


STRAIN RANGE 

15. 000 
S.R000 

1.0000 

0. 10000E+00 


STRESS* TO«RUPTllRE CURVE 


HOURS 

0,]10000E + 00 
0,60000 

6,0000 

100,00 

1000.0 


STRESS (KSI) 

6.0000 

4.0000 

2.5000 

1.5000 

1 .0000 



WtSuLTS FQN Each CYCLt 


U1 


NCYC 

tmi m 

tmax 

EtQvT 

N/Ntl 

l 

0,ET80lE-01 

0.28093E-01 

2. *630 

0,22575t-02 

2 

0,27b03E-Ol 

0.28187E-01 

2.9169 

0.231 Jlfc-02 

: * 

0.2740«E-01 

0.282806-01 

2.9517 

0.23708E-02 

a 

0,’2720bE-0l 

0.28374t-01 

2.9879 

O.24308E-02 

5 

0,27007E-01 

0.26467E-01 

3.0241 

0.24932E-02 

6 

0 ,26809E-0 1 

0. 285616-01 

3.0617 

0.25S81E-02 

7 

0,2661oE-01 

0.28654E-01 

3.1003 

0,2b25bE-02 

8 

0,2641 26-0 1 

0,287476-01 

3.1400 

0,2b959t-02 

9 

0,2621 3E-01 

0.2B841E-01 

3.1807 

0.27690E-02 

10 

0. 2601 SE-0 1 

0.289.346-01 

3.2224 

0.28451E-02 

It 

0^25616E-01 

0.2902AE-01 

3.2653 

0.29243E-02 

12 

0,25blsE-01 

0.29121E-01 

3.3093 

0.30067E-02 

13 

0,254196-01 

0.29215E-01 

3.3544 

0.30926E-02 

14 

0,2522lE«01 

0.29308E-01 

3.4007 

0.31820t-02 

15 

o,25oaaE-oi 

0.294016-01 

3.4482 

0 .327 516-02 

16 

0,24B24E-01 

0.29495E-01 

3.4970 

0. 337216-02 

17 

0,2«625E-01 

0 . 29588E-0 1 

3.5470 

0.34731E-02 

16 

0,2442bE-01 

0.29682E-01 

3.5984 

0.35784E-02 

19 

0,2422flE-0l 

0.297756-01 

3.6511 

0 , 36882E-02 

20 

Q,24029E-01 

0,29B69E-01 

3.7051 

0 .380266-02 

21 

0,238J]E-01 

0.29962E-01 

3.7606 

0 , 3921 9t-02 

22 

0,236326-01 

0.300556-01 

3.8175 

0,404b2t-02 

23 

(J.2343«E-0l 

0.301496-01 

3.8759 

0.41759E-02 

24 

0,2323sE»01 

0.30242E-01 

3,9358 

0,431126-02 

25 

0,23037E»01 

0.30336E-01 

3.9973 

0.44523E-02 

26 

0,22a3sE-01 

0, 30429E-01 

4,0604 

0 ,459966-02 

27 

"0,226«6E»6~1 

0.30522E-01 

4.1251 

0 ,475336-02 

26 

0,224«iE-01 

0.30616E-01 

4.1915 

0,491376-02 

29 

0,2224jE-0l 

0.307096-01 

4.259b 

0,506126-02 

30 

0,220446-01 

0.30803E-01 

4.3296 

0.52560E-02 

31 

0,2t84bE-01 

0.30896E-01 

4,4013 

0.54387E-02 

32 

0,2164TE-01 

0.30990E-01 

4.4749 

0 . 56295E-02 

33 

6,2U49EVoi 

0.310B3E-01 

4.5505 

0 ,582886-02 

34 

0,21250E-01 

0.31176E-01 

4.6280 

0,6037 1 E-02 

35 

0,2 1 05l E«0 1 

0.31270E-01 

4.7076 

0,625486-02 

36 

0,20853Ef0l 

0,313636-01 

4.7892 

0,64824E-02 

37 

0,2065flE«01 

0.31457E-01 

4.8731 

0,672046-02 

36 

0«g0456E«Ql 

0,315506-01 

4.9591 

0.69693E-02 

34 

0,2025?E«0l 

0.31644E-01 

5.0474 

0 , 72297E-02 

40 

0,20059E-01 

0.31737E-01 

5.1380 

0,750216-02 

41 

Q l 19860E-01 

0.318306-01 

5.2311 

0.77872E-02 

42 

0,1966aE-01 

0.31924E-01 

5.326b 

0,808576-02 

43 

0,14«6jE-0t 

0.32017E-01 

5.4247 

0,839826-02 

44 

0,1926SE«01 

0,3211 IE-01 

5.5254 

0.87255E-02 

45 

6,19066E-dl 

0, 32204E-0 1 

5.6288 

0.90683E-02 

46 

0 , 18868E-0 1 

0.32298E-01 

5.7350 

0,942756-02 

47 

0,18o69E-01 

0.32391E-01 

5,8441 

0 , 9804 1 E-02 

48 

0,1847 1 E w 0 1 

0.32484E-01 

5.9561 

0.10197E-01 

49 

0, 1827jE«01 

0,325786-01 

6.0712 

0,1 06066-01 

5Q 

0»180T4E**01 

0,3267 16-01 

6,1895 

0,1 10396-01 


STK6S3 

r/ro 

TOTAL USAGE 

NF UYCLtS 

0.0000 

o.oooo 

0.225756-02 

442 

0,0000 

0,0000 

0.457066-02 

432 

0.0000 

o.oooo 

0,694146-02 

421 

0.0000 

O.oooo 

0.93722E-02 

411 

0.0000 

O.oooo 

0,118656-01 

401 

0,0000 

0,0000 

0.144246-01 

391 

0,0000 

0.0000 

O.17049E-01 

381 

0.0000 

0,0000 

0.197456-01 

371 

0.0000 

0,0000 

0 , 225 1 4E-0 1 

362 

0.0000 

0,0000 

0.253596-01 

352 

0.0000 

0.0000 

0,282836-01 

343 

0.0000 

0,0000 

0.312906-01 

314 

0,0000 

0.0000 

0 , 34383E-0 1 

32b 

0.0000 

0,0000 

0.37565E-01 

316 

0,0000 

0.0000 

0.40840E-01 

307 

0,0000 

0.0000 

0.442126-01 

299 

0.0000 

0,0000 

0.47685E-01 

291 

0,0000 

0,0000 

0.51263E-01 

283 

0,0000 

0.0000 

0.54952E-01 

275 

0,0000 

0.0000 

0,587546-01 

267 

0.0000 

0.0000 

0 .62676E-0 1 

259 

0,0000 

0,0000 

0.66722E-01 

252 

0,0000 

0,0000 

0.70898E-01 

245 

6,0000 

o.oooo 

0,'75209E-01 

236 

0.0000 

0.0000 

0 , 79662E-0 1 

231 

0.0000 

0.0000 

0,842616-01 

225 

6,0000 

O.oooo 

0.890156-01 

216 

0,0000 

o.oooo 

0.939286-01 

212 

0,0000 

0.0000 

0.99009E-01 

206 

0,0000 

0,0000 

0,10427 

200 

0.0000 

0,0000 

0,10970 

194 

0^0000 

0.0000 

0.11533 

189 

0,0000 

O.oooo 

0.12116 

183 

0,0000 

0.0000 

0.12720 

176 

0.0000 

0,0000 

0.13345 

173 

0,0000 

0,0000 

0,13994 

166 

0,0000 

0,0000 

0,14666 

163 

0,0000 

0,0000 

0,15363 

159 

0,0000 

6,0000 

0,16086 

155 

0,0000 

0,0000 

0.16836 

150 

0,0000 

0,0000 

0.17615 

146 

6,0000 

O,o000 

0^18423 

142 

0,0000 

O,0000 

0,19263 

139 

OjOOOO 

0.0000 

0.20135 

135 

0,0000 

6.6000 

o;21042 

132 

0,0000 

0,0000 

0,21985 

126 

0,0000 

0,0000 

0.22965 

125 

6,0000 

0,0000 

0.23985 

122 

0,0000 

0,0000 

0,25046 

119 

0.0000 

n.0000 

0.26130 

116 



cr> 


51 

0 t lT875E-01 

0.32765E-01 

6,3110 

0,1 14926-01 

0,0000 

0,0000 

0,27299 

52 

0^lTfa7bE-01 

0.32858E-01 

6.4359 

0,119666-01 

0,0000 

0,0000 

0.28496 

53 

0,1TU78E-01 

0.32952E-01 

6.5643 

0.12464E-01 

o.oooo 

0,0000 

0.29742 

50 

0,lT27qE»0i 

0.33045E-01 

6.6962 

0, 12986E-01 

0,0000 

0,0000 

0.31041 

55 

oJlTOBlE-Ol 

0 . 33 138E-01 

6.8319 

0,1 353SE-0 1 

0.0000 

C.0000 

0.32394 

56 

0,1688aE-01 

0.33232E-0! 

6.9714 

0.141126-01 

0.0000 

0,0000 

0,33606 

57 

o Jl668oE-01 

0.33325E-01 

7.1108 

0,147186-01 

o.oobo 

0,0000 

0.35277 

58 

0,16485E-01 

0.33419E-01 

7.2624 

0 . 1 5355t-0 1 

0,0000 

n.oooo 

0.36813 

59 

0 jl6287E-01 

0.335J2E-01 

7.0102 

0.16025E-01 

0,0000 

^ * oooo 

0.58415 

60 

O.lf'OSsE-Ol 

Cr.33606h-0l 

03 

0.167306-01 

0,0000 

0,0000 

0*40066 

61 

o,i589oE-oi 

O.33699E-01 

7.7311 

0,174716-01 

0.0000 

0,0000 

0.41835 

62 

0 ,1569 1 E-01 

0.33792E-01 

7.8966 

0, 182526-01 

0.0000 

0,0000 

0,43661 

63 

0,15493E-0i 

0.33886E-01 

8.0670 

0.190756-01 

0,0000 

0,0000 

0 , 45566 

60 

0,15294E-01 

0.33979E-01 

8.2424 

0,199416-01 

0,0000 

0,0000 

0.47562 

65 

0*1509bE-01 

0.34073E-01 

6.4232 

0.208546-01 

0.0000 

0,0000 

0.49646 

66 

0 t 14897E*01 

0.34166E-01 

8.6094 

0,218l7E-01 

0,0000 

0,0000 

0,51829 

67 

0,14699E-01 

0.34260E-01 

8.8013 

0,228336-01 

0,0000 

0,0000 

0,54113 

66 

0,t4500E-01 

0 , 34353E-0 1 

8.9991 

0.23905E-01 

0,0000 

0,0000 

0.56503 

69 

b,i43b{E-oi 

0.34446E-01 

9.2031 

0,2S037E-01 

0,0000 

o.oooo 

0,59007 

70 

0.14103E-01 

0.34540E-01 

9.4135 

0,262336-01 

0,0000 

0,0000 

0.61630 

71 

Q,13904E-01 

0.34633E-01 

9.6306 

0,274976-01 

0.0000 

0.0000 

0.64380 

72 

0,13706E-01 

0.34727E-01 

9.8546 

0,288346-01 

0,0000 

o,oooO 

0.67263 

73 

0.13507E-01 

0.34820E-01 

10.086 

0,302496-01 

0,0000 

0,0000 

0,70288 

74 

75 

0* J3309E"0l 
o,i3ii6£-oi 

0.34913E-01 

6»350b7E“01 

10.325 

10.572 

0.31746E-01 
0 .333326-01 

0,0000 

o.oobo 

0,0000 

b.oboo 

0,73463 
0. 76796 

76 

0,12912E-01 

0.35100E-01 

10.827 

0,350146-01 

0,0000 

0,0000 

0.80297 

77 

0,12713E-01 

0.35194E-01 

11.090 

0.36797E-01 

0,0000 

0,0000 

0,83977 

76 

0.12515E-01 

0.35287E-01 

11.363 

0,386896-01 

0,0000 

0,0000 

0,67846 

79 

0.12316E-01 

0.35381E-01 

11,645 

0,406996-01 

0,0000 

0,0000 

0,91916 


C v 12116E-01 

0,354746-01 

11,937 

0,428346-01 

0,Q000 

0.0000 

0,96199 

61 

0.11919^-01 

0.35567E-01 

12.240 

0,451056-01 

0,0000 

o.oooo 

1.0071 


114 

111 

109 

107 

104 

102 

100 

99 

97 

9b 

99 

92 

91 

90 

89 

86 

87 

66 

85 

64 

83 

83 

62 

62 

61 

81 

61 

61 

60 

80 

60 



CONSTANTS INPUT 


Qa 


0.22*00 


DELTA 

■ 

0.37500E-03 

INCHES 

E 1 A VG 


W4250 

PERCENT 

H 


0.14000E-01 

INCHES 

L 


O.OOOOOE-Ol 

INCHES 

W 


o.asoonE-oi 

INCHES 

P 


3000.6 

LBS/INCH 

T/CYC 


4.0000 

MINUTES 

SF 

a 

o.aoooo 



fatigue cuhve 


CYCLE? 

jo.ooo 

100,00 

000*0 

0. 5Q000E*0b 


STRAIN RANGE 

IB.’OOO 
5* ROOD 

1.0000 

0, 10000E+00 


STRESS-TO-RUPTURE CURVE 


HOURS 

0^1 OOOOE+OO 
0,80000 

8,0000 
100. Qo 

1000.0 


STRESS (KSI) 

6.0000 
a. oooo 

2.5000 

1.5000 

1.0000 



00 


Rt SUL T S FOR EACH CYCLE 


NCYC 

THIN 

TM*x 

FfcQvT 

N/fcO 

STRESS 

T/TO 

TOTAL usage 

Nf CYCLES 

t 

0 , '£ 78016-01 

0.28093E-01 

2.8830 

0.225756-0? 

3.0323 

0.214556-01 

0.23713E-01 

42 

a 

O^ToOjE-Ol 

0.28187E-01 

2.9169 

0.2313 1 £-02 

3.0761 

0 , 230166-01 

0,49042E-01 

39 

3 

0,274046-01 

O.2828OE-01 

2.9517 

0.23708E-02 

3.jaofl 

0.247036-01 

0.76116E-O1 

3 1 

Q 

0,2720bE-01 

0.2857 4E-0 1 

2.9874 

0 . 24308E-0? 

3.1665 

0.265276-01 

0.10507 

34 

5 

0,270076-01 

0.28467E-01 

3.0241 

0.249326-02 

3 . ? 1 33 

0.285016-01 

0,13607 

32 

6 

0,26B09E-01 

0 .2856 IE-01 

3 . Ob 1 7 

0.255816-02 

3.2610 

0.306376-01 

0.16926 

3i 

7 

0 ,266 1 OE-0 1 

0.28654E-01 

3.1003 

0.2625(«E-02 

3.3099 

0.329526-01 

0.20484 

29 

ft 

0,264lpE-01 

0.28747E-01 

3.1400 

0.2b959E-02 

3.3598 

0,354616-01 

0,24300 

27 

9 

0,2621 3E-01 

0.23841E-01 

3.1807 

0.276906-02 

3.4109 

0,381816-01 

0.28395 

26 

10 

0,260 1 SE -0 1 

0 . 28934E-0 1 

3.2224 

0.284516-02 

3.4632 

0,41 1346-01 

0.32793 

25 

ii 

0,258 1 feE-0 1 

0.2902BE-01 

3.2653 

0.292436-02 

3*5188 

0,443416-01 

0,37519 

24 

12 

0,25b 1 dE-0 1 

0.291 21E-01 

3.3093 

0.300676-02 

3.5714 

0.478256-01 

0.42602 

23 

13 

0,254 196-01 

0.2921SE-0t 

3.3544 

0.309266-02 

3.6274 

0. 516136-01 

0,48073 

22 

14 

0 ,2522 1 E-0 1 

0.29308E-01 

3.4007 

0.318206-02 

3.6647 

0.55734E-01 

0,53964 

21 

15 

0,25022E-0l 

0.29401E-01 

3.4482 

0.327516-02 

3.7434 

0 , 602216-01 

0,60314 

21 

16 

0,20824E-01 

0.29495E-01 

3.4970 

0,337216-0? 

3.8035 

0 . 651106-01 

0.67162 

20 

17 

0,246256-01 

0.29588E-01 

3.5470 

0 • 34731 6-02 

3.9651 

0,704406-01 

0,74553 

20 

1 ft 

0,2442bE-01 

0.29682E-01 

3,5984 

0.357846-02 

3.9282 

0,762546-01 

0,82537 

20 

19 

0,24228E-01 

0 . 29775 Er 01 

3.6511 

0,368826-02 

3.99?B 

0 , 826026-01 

0.91166 

20 

20 

0.2402<?E-01 

0.29869E-01 

3.7051 

0.380266-02 

4.0591 

0.898386-01 

1,0053 

19 



APPENDIX C 
FORTRAN PROGRAM FOR 

FATIGUE AND CREEP RUPTURE DAMAGE EVALUATION 
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ocio onci ooorioo nnnnnn nnricionoooooor* 


PROGRAM CYCLES 

REAL NF f l.,No 

DIMENSION CyCLE(20) ,STRAIN(20) ,STG(20) , TIME (20) 

C« EXPONENT in STRESS-STRAIN lA*. 

DELTA. OEEOrMATION P EH CYCLE 
E 1 A VG . AVERAGE HOUP STRAIN (PERCENT) 

H. HAtF TME HEIGHT OF LIGAMtNT 
L. LENGTH OF ligament 
W. WIDTH OF RIR 
P. PRESSURE (LBS/IN) 

TPC » TIME PER CYCLE (MlNUTFS) 

SF. shape factor 

TU ELIMINATE CREEP EFFFCTS,' ENTER 0 FOR P, TPC AM) SF 
RE AO (15,*) Q»OElTA,’tlAVG,H,L.'w,P,TPC,SF 

input fatigue curve 

nfat.nukber of data points on fatigue curve to be input (max of 20 ) 
CTCLE.CYCLEs for each data pt 
strain. STRAIN RANGE FOR EACH PT 

READ ( 15, *) NEAT 

READ ( 15, *) (CYClE(I),I«i,NFAT) 

READ ( 15, * ) (STRAIN(I) ,I»1,NFAT) 

INPUT ST RE $S*T 0-RUPTURE CURVE 

NRUP.NUMBER OF DATA POINTS TO BE INPUT (MAX OF 20) 

SIG. STRESS FOR EACH DATA PT (KSI) 

time. time for each data pt (hours) 

READdS,*) nRUP 

READ(15,*) (TIME(n,I«l,NRUP) 

REA0(15,*) (SIGCD.’IM.NRUP) 

PRINT SUMMARY OF INPUT 

WRITE (16, 500) Q,0ELTA,E1AVG,H,L,W,P,TPC # SF 
DO 5 I« 1 ,NFaT 

5 WRITE(16,600) CYCLEII) ,'STRAINCI) 

WRITE(16,620) 

DO 7 I • 1 » NRuP 

7 WRITE (lb, 6 O 0 ) TIME ( I ) , SIG ( I ) 

calculate cycles to failure 

WRI TE ( 1 6 , 7 00) 

USEAGE»0»0 
NCVC»0 
C»(0-l)/0 
TPC.TPC/bO 
10 NC YC»NC YC* 1 

TMIN»(2*H.(l*W)-NCYC*0FLTA.W)/CL*W) 
TMAX«(2*H*(L*H)**2*NCTC*0ELTA*L*W) /CL*W)**2 
T»TMAX/TMIN 

ElMIN«ElAVG.C*(T-I)/n**C-n 
E 1 *E 1 MIN 
E2»E 1 A VG 

EEOVT.SORT ( (E l-E?) **2* (2*E 1 *E2) **2* (2*E2*t 1 ) *.2) *SOKT (2 , ) /5 

00 20 I *2 , Np A T 

IF(EE0VT,LE.STRAIN(I-1 ) .and. EEOVT.GE.STRAIN(D) goto 30 
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nno non 


20 CONTINUE 

WRITE ( 16*300) EEOVT 

30 No 2 exP(ALOG(EEQVT/ 9 TR* IN (I-l))/AL 0 G (STRAIN (I) /STRAIN ( 1 - 1 ))* 

2A10G (CYCLE (I) /CYCLE (I-1))*AL0G (CYCLE (I-l>)> 
calculate TIME TO ruptuhe 


TD»1 

IF(P.EO.O.O) GO TO 60 
PM»SF*P* (L*W/2)"*2/24 
Sb6*PH/TM1N**2 


S-S/1000 
00 40 I"2, NRUP 

IF(S.LE.3IG(I»1) .AND, S.GE.SIG(I)) GO TO 50 
40 CONTINUE 

IP(S.LT.SIG(NRUP)) GO TO 45 

T0"EXP(AL0G(TIME(2))-AL0G(9IG(2)/S)/AL0G(3IG(2)/SIG(1) 

2 AL0G(TlME(2)/TIME(im 


)* 


GO TO 60 

45 T 0 »EXP(AL 0 G(TIHE(NRUP- 1 ) )«ALOG (SlG(NRUP-n/S) / 

2 ALOG(SIG(NrUP*1)/SIG(NRUP)) *alog (time (nrup*D /time CNRUP))) 


50 TD«EXP (ALOG (S/S1G(I"1 ) ) / ALOG(SlGd) /SIG (1*1 ) ) * 
2 ALOG(TIME(l)/TIME(I-J))+ALOG(TIME(I-l)n 


compute cumulative DAMAGE 


60 useage«useage ♦ l./NO ♦ TPC/TO 

NF»(l«USEAGE)/(l./NO ♦ TPC/TD) ♦ N(JyC 
WRITE (16,600) NCYC,TMIN,TMAX,EE0VT,1. 
IF(USEAGE.LT.l.O) GO TO 10 
STOP 


/N0,3,TPC/ID, 


USEAGE, IFIX (NF) 


FORMATS 

TOO FORMAT(1 h O, "THE STRAIN RANGE OF ",Gl2.5, 

1 » IS NOT ON THE FATIGUE CURVE SUPPLIED") 

350 FORMATdHO, "THE STRESS OF ",Gl2 ( 3i 

1 • IS NOT ON THE 3TRESS-T0-RUPTURE CURVE SUPPLIED ) 

500 FORMAT d H l , 12X /"CONSTANTS INPUT"// 1 1 X , " 0" ,G12,5//11X, 

2 "OELTA "**612,5," INCHES"//11X, 

PERCENT*//11X, 

INCHES*//! IX, 

1NCHES*//HX, 

INCHES"//1 IX, 

lbs/inch"//i IX, 


3 

4 

5 

6 
7 
6 
6 


•ElAVG "",012. 5,.' 

» H «»,G12.5*" 

* L ■'',612. 5„* 

» W »*,Gl2.5,," 

P "»,G12,5;* 

"T/CVC "*,G12,5,» Ml NUTF3" / / 11 X , 

„ " SF "",G12.5////18X, 

9 "Fatigue cuRve b //i4x, "cycles", 7x, "strain range"/) 
600 FORHATdH , b* , 2 ( «X ,G1 2 , 5) ) „„ „ „ 

620 FORMA T (1 HO/// 13 X," 9 TRESS"T 0 “RUPTURE CuRvE"//13X, 

2" HOURS", 6X, "STRESS (NSI)"/) 

700 FORMAT d H l , 46X , "RESULTS FOR EACH CYCLE"// 

2 

3 

800 FORHATdH ,I7,7(3X,G12.5) ,18) 

END 


•0RHATdHl,46X, "RESULTS TOR EACH 0 YCLt / , 

5X,"NCYC*,6X, "THIN", 11X,*TMAX",11X, "EEOVT", 10X, N^D , 
10xV"3TRESS",10X, M T/T0",9X, "TOTAL USAGE ", TX , "NF CYCLES /) 
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APPENDIX D 


SYMBOLS 

A - area 

"A - equation arbitrary positive scalar 

B - material constant in Norton's creep law 

D - constant in Ramberg-Osgood stress strain law 

e - elastic strain 

E - modulus of elasticity 

JF (r ) - Factor in Figure 10 

F(r) - Factors in Table I 

2H - thickness of ligament 

k - /I - s z ; also thermal conductivity 

K - curvature 

l - width of ligament in hoop direction 

m - generalized bending stress variable 

M - bending moment in ligament 

M 0 - yield bending moment 

n - generalized hoop stress variable 

N - hoop force in ligament; also number of cycles 
N q - yield hoop force 

p - pressure 

q - exponent in Ramberg-Osgood stress strain law 

r - exponent in Norton's creep law 

s - generalized shear stress variable 

S - shear force in ligament 

Sy - average yield stress in tension 

S - ligament yield strength for minimum a(T.j - T ) 
y min 

S - ligament yield strength for maximum a(T- - T ) 
y max 0 
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SYMBOLS - continued 


cr 


rmn 


max 


'ss 


T 

T, 


w 


time 

critical ligament thickness 

thinning after N cycles 

minimum ligament thickness 

maximum ligament thickness 

time to reach steady state distribution 

temperature 

average temperature of ligament 
average temperature of closeout wall 
width of rib 


a 

Y 

6 

61 

<$2 

63 




El 

El 


£l avg " 
El min 


mi 

cr 


coefficient of thermal expansion; also stress ratio a 2 /ai 

shear strain 

deflection per cycle 

deflection due to moment 

deflection due to shear 

creep induced deflection 

pressure difference between coolant pressure and combustion 
gas pressure 

inelastic strain range in hoop direction due to differential 
thermal expansion 

correction to plastic strain range in hoop direction due to 
thermally induced bending 

strain 

hoop strain 

average hoop strain in ligament 

hoop strain in minimum ligament section 

axial strain in minimum ligament section 
critical effective strain 



SYMBOLS - conti nued 


£ min " effective in minimum ligament section 

0 - generalized bending strain variable 

X - generalized hoop strain variable 

v - Poisson’s ratio 

a - stress 

01 - hoop stress 

a 2 - axial stress 


Aa 


a in 


a 


ss 


<f> 


creep relaxation stress 

initial pressure induced bending stress 

peak stress 

steady state stress 

dimensionless time 

generalized shear strain variable 

dot above symbol denotes rate 
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